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, Abstract : In [Sj we have denned a viscosity solution for the gradient 

flow of the exterior Bernoulli free boundary problem. We prove here that 
the associated energy is non decreasing along the flow. This justifies the 
"gradient flow" approach for such kind of problem. The proof relies on the 
construction of a discrete gradient flow in the flavour of Almgren, Taylor 
and Wang |2j and on proving it converges to the viscosity solution. 

> 

in 

■ 1 Introduction 
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. free boundary problem initiated in |S]. The exterior Bernoulli free boundary 

problem amounts to minimize the capacity of a set under volume constraints. 
Using a Lagrange multiplier A > 0, this problem can be recasted into the 
minimization with respect to the set f2 of the functional 
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In this paper we continue our investigation of a gradient flow for the Bernoulli 



£\(n) = cap 5 (ft) + A|0| 

■ where ca,p s (Q) denotes the capacity of the set Q with respect to some fixed 

set S and |f2| denotes the volume of O. The set £1 is constrained to satisfy 
the inclusion S CC 0. Notice that there is a "competition" between the 
two terms in the minimization: the capacity is nondecreasing with respect 
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to inclusion whereas the volume is nondecreasing. 



Such a problem has quite a long history and we refer to the survey paper 
|12| for references and interpretations in Physics. Our study is motivated by 
several papers in numerical analysis where discrete gradient flows are built 
via a level-set approach in order to solve free boundary and shape optimiza- 
tion problems: see ^ and the references therein for the recent advances in 
this area. In this framework, the exterior Bernoulli free boundary problem 
appears as a model problem in order to better understand this numerical 
approach. In this work, we prove that the energy £\ is non increasing along 
the generalized flow we built in jSJ. This question is certainly essential to 
better explain the numerical schemes of pQ. This also fully justifies the ter- 
minology of "gradient flow" for the generalized solutions. 

Let us now go further into the description of the gradient flow for £ := £\ 
(we work here in the case A = 1 for simplicity of notations). The energy 
£ being defined on sets, a gradient flow for £ is a family of sets (Q(t)) t >o 
evolving with a normal velocity which "decreases instantaneously the most 
the energy" . For the Bernoulli problem, the corresponding evolution law is 
given by: 

V ttX = h(x,Q(t)) := -l + h(x,n(t)) for all t > 0, x G dfi(i) . (1) 

In the above equation, Vj )IE is the normal velocity of the set f2(i) at the point 
x at time t and h(x, f2) is a non local term of Hele-Shaw type given, for any 
set O with smooth boundary, by 



where u : £1 — > R is the capacity potential of ft, with respect to S, i.e., the 
solution of the following partial differential equation 

( -Ait = in 9\S, 



The set S is a fixed source and we always assume above that S is smooth 
and S CC fi(t). Let us underline that h(x,£l) is well defined as soon as f2 
has a "smooth" (say for instance C 2 ) boundary and that S CC O. 

The reason why a smooth solution (O(t)) of the geometric equation Q 
can be considered as a gradient flow of the energy 



h(x,Q) = \Vu(x)\ 
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(2) 




on dS, 
on dU. 



(3) 



£(n) = \n\ + cap s (f2) 



(4) 
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is the following: from Hadamard formula we have 




Hence the choice of V tjX = h(x,Q(t)) in (JJ) appears to be the one which 
decreases the most the energy £. In order to minimize the energy £, it is 
therefore very natural to follow the gradient flow (^Q). This is precisely what 
is done numerically in 1 . 

In general the geometric flow Q does not have classical solutions. In 
order to define the flow after the onset of singularities, we have introduced in 
[S] a notion of generalized (viscosity) solution and investigated its existence 
as well as its uniqueness. In order to prove that the energy is non increas- 
ing along the generalized flow, we face a main difficulty: energy estimates 
are hard to derive from the notion of viscosity solutions. Indeed this latter 
notion is defined through a comparison principle, which has very little to 
do with the energy associated to the flow. To the best of our knowledge, 
such a question has only be settled for the mean curvature motion (MCM 
in short), which corresponds to the gradient flow of the perimeter. There 
are two proofs of the fact that the perimeter of the viscosity solution to the 
mean curvature flow decreases: the first one is due to Evans and Spruck 
in their seminal papers |10[ 111): it is based on a regularized version of the 
level set formulation for the flow and is probably specific to local evolution 
equations. The other proof is due to Chambolle 01 • Its starting point is 
the fondamental construction of Almgrem, Taylor and Wang 2 who built 
generalized solutions of the MCM in a variational way as limits of "discrete 
gradient flow" for the perimeter (the so-called minimizing movements. See 
also Ambrosio 5 ]). The key argument of Chambolle's paper j^j is that Alm- 
gren, Taylor and Wang's generalized solutions coincide with the viscosity 
solutions, at least for a large class of initial sets. Hence the energy estimate 
available from (2j — which allows to compare the energy of the evolving set 
with the energy of the initial position — can also be applied to the viscosity 
solution. Since the viscosity solution enjoys a semi-group property, one can 
conclude that the energy is decreasing along the flow. 

For proving that the energy £ is decreasing along our viscosity solutions 
of we borrow several ideas from Almgren, Taylor and Wang j2j and 
Chambolle [S]. As in for the MCM, we start with a construction of dis- 
crete gradient flow (O^) for the energy £: namely is obtained from 
as a minimizer of a functional Jh(Q^, •) which is equal to £ plus a penalizing 
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term. The penalizing term — which depends on the time-step h — prevents 
the minimizing set Q^+i fr° m being too far from Cl^. Then, as in Chambolle 
[S] , we prove that the limits of these discrete gradient flows converge to the 
viscosity solution of our equation (JJ) as the time-step h goes to 0. In 
this convergence is proved by using the convexity of the equivalent of our 
functional ,7^(0^, •) for the MCM. We use instead here directly a weak form 
of the Euler equation for minimizers of J^QV;,-) as described by Alt and 
Caffarelli j3j for the Bernoulli problem. We then conclude that the energy 
of the flow is non increasing. 

The paper is organized in the following way. In Section [2 we recall the 
construction of jS] for the viscosity solutions of (JJ) . Section |3] is devoted 
to suitable generalizations of the capacity and capacity potential needed for 
our estimates. In Section @] we introduce the functional Jh and build the 
discrete motions, the limits of which are discussed in section The fact 
that the energy is decreasing along the flow is finally proved in Sectional 

Aknowledgement : We wish to thank Luis Caffarelli, Antonin Cham- 
bolle and Marc Dambrine for fruitful discussions. The authors are partially 
supported by the ACI grant JC 1041 "Mouvements d'interface avec termes 
non-locaux" from the French Ministry of Research. 

2 Definitions and notations for the generalized flow 

Let us first fix some basic notations: if A, B are subsets of R N , then A CC B 
means that the closure A of A is a compact subset which satisfies A C int(5), 
where int(B) is the interior of B. We set 



Throughout the paper | • | denotes the euclidean norm (of W N or R^ 1 , 
depending on the context) and B(x,R) denotes the open ball centered at x 
and of radius R. If E is a measurable subset of W N , we also denote by \E\ 
the Lebesgue measure of E. If K is a subset of M. N and x € M. N , then dx(x) 
denotes the usual distance from x to K: (1k(x) = inf yg A' \y — x\. The signed 
distance d s K to K is defined by 



V = {K CCR N : S CCK} . 
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where dK = K\int(K) is the boundary of K. Let Vt be an open bounded 
subset of R N . We denote by C£°(fi) the set of smooth functions with com- 
pact support in Q, and by Hq(Q) its closure for the H 1 norm. 

Here and throughout the paper, we assume that 

S is the closure of an open, nonempty, bounded subset of R N , , 
with a C 2 boundary. 

The generalized solution of the front propagation problem (JJ) is defined 
though their graph: if (Q(t))t>o is the familly of evolving sets, then its graph 
is the subset of R + x R N defined by 

K. = {(t,x) £R + xR N : x G Sl(t)} . 

We denote by (t, x) an element of such a set, where t G R + denotes the time 
and x G R N denotes the space. We set 

K{t) = {x e R N | (t,x) € K.} . 

The closure of the set K in M^ 4 " 1 is denoted by K. The closure of the 
complementary of K is denoted K: 

t = (R+ x R N ) \/C 

and we set 

= {x G | (t,x) G £} . 
We use here repetitively the terminology of (HI El IH1 : 

• A tube K, is a subset of M+ x R N , such that ^ n ([0,i] x E^) is a 
compact subset of M Ar+1 for any i > 0. 

• A tube /C is Ze/t lower semi- continuous if 

Vi > 0, Vx 6 /C(t), if t n — > 3x„, G K.{t n ) such that x n — > x . 

• If s = 1, 2 or (1, 1), then a C s tube K, is a tube whose boundary <9/C 
has at least C s regularity. 

• A regular tube fC r is a tube with a non empty interior and whose 
boundary has at least C 1 regularity, such that at any point (t, x) G /C r 
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the outward normal (z/f, u x ) to fC r at (t, z) satisfies v x 7^ 0. In this 
case, its normal velocity Vu r x \ a t the point (t, x) G dfC r is defined by 

T/^ r = — 
(*.*) |i/ a |' 

where {yti v x) ls the outward normal to /C r at (£, z). 

• AC 1 regular tube )C r is externally tangent to a tube /C at (i, x) G XI if 

K. <Z K. r and (i, z) £ <9/C r . 
It is internally tangent to /C at (i, z) G /C if 

/C r C K and (i, z) G dtC r . 

• We say that a sequence of C ' tubes (fC n ) converges to some C ' tube 
/C in i/ie C 1,b sense if (/C n ) converges to K and (dlC n ) converges to 
dlC for the Hausdorff distance, and if there is an open neighborhood 
O of dK such that, if a\ (respectively a\ ) is the signed distance to 
K, (respectively to K n ), then (df^ ) and (Vdj^ ) converge uniformly to 
<ijc and -Dd^ on and \\D 2 d s jQ n are uniformly bounded on O. 

We are now ready to define the generalized solutions of (J2): 

Definition 2.1 Let fC be a tube and Kq G T> be an initial set. 

1. K is a viscosity subsolution to the front propagation problem (QJ) if JC 
is left lower semi- continuous and IC(t) G T> for any t, and if, for any 
C 2 regular tube JC r externally tangent to K, at some point (t,x), with 
IC r (t) G T> and t > 0, we have 

V^ x) <h{x,Kr{t)) 

where is the normal velocity of K r at (t,x). 

We say that fC is a subsolution to the front propagation problem with 
initial position Kq if JC is a subsolution and i//C(0) C Kq. 

2. tC is a viscosity supersolution to the front propagation problem if fC 
left lower semi- continuous, and IC(t) C V for any t, and if, for any 
C 2 regular tube fC r internally tangent to fC at some point (t,x), with 
/C r (i) G V and t > 0, we have 

Vfc x) > h {x,1C r {t)) . 

We say that K, is a supersolution to the front propagation problem with 
initial position Kq if fC is a supersolution and if K(0) C M. n \Kq. 
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3. Finally, we say that a tube /C is a viscosity solution to the front prop- 
agation problem (with initial position Kq) if K is a sub- and a super- 
solution to the front propagation problem (with initial position Kq). 

In [Sj we have proved that for any initial position there is a maximal 
solution, with a closed graph, which contains any subsolution of the problem, 
as well as a minimal solution, which has an open graph, and is contained in 
any supersolution of the problem. 

3 Capacity and capacity potential 

Let S be as in ©. For an open bounded subset Q of M. N such that S CC 
the capacity of Q with respect to S is defined by 



Since S is a fixed set in what follows, we will write cap(J7) instead of cap 5 (J7). 

Obviously cap(O) is non increasing with respect to the set £1 (for inclu- 
sion). For a general reference on the subject, see for instance |14j . 



Remark 3.1 (Classical capacity potential) If f2 is any bounded open subset 
of R N , then 



and the infimum is achieved for a unique u £ Hq(Q), called the capacity 
potential of ft, with respect to S, such that u = 1 on S, u is harmonic in 
Q\S and \{u > 0}\O| = (namely, u = a.e. in R N \Q). If O has a 
C ' boundary, then it is known that the infimum is achieved by a function 
u G C 2 (Q\S) n C x (r2\5) which is a classical solution to ©• 

For any set E (not necessarily open) such that S CC E, we define a 
generalized capacity by 

cap(E) = sup{cap(0) | E CC £1, 0, open and bounded} . 

With this definition, cap(-E) is non increasing with respect to the set E. 
Notice that this notion of capacity does not take into account "thin closed 
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sets" in the sense that, if F = E, then cap(E) = cap(-F) even when ^ 
0. By construction, if E is open, then we have 

cap(-E) < cap(E') 

but equality does not hold in general. Nevertheless, there is equality if the 
boundary of the set is regular enough: 

Lemma 3.1 If Q, is an open bounded subset ofM> N , with S CC SI and with 
a C ' boundary, we have cap(f2) = cap(f2). 

Proof of Lemma 13.11 We have to prove that cap(f2) > cap(f2). It is 
enough to show that, if 

n n = { y em N : d n (y) < 1/n} , 

then cap(f2 n ) — > cap(f2) as n — > +oo. Indeed, for n large enough, f2 n 
has also a C 1,1 boundary. Then from classical regularity arguments, the 
harmonic potential u n to Q n converges to the capacity potential u of O for 
the C 1,a norm, where a S (0, 1). Whence the result. 

QED 

Lemma 3.2 Let E n be a bounded sequence of subsets of~R N , for which there 
exists some r > with S r C E n for any n, where 

S r = {y e M N : d s (y) < r} . (7) 

Let us denote by K the Kuratowski upper limit of the (E n ), namely 

K = jx G 1^ : lirninf d En (x) = o} . 

Then 

liminf cap(£' n ) > cap(K) . 

n 

Proof of LemmalO Let tt be any open bounded set such that K CC fl. 
Since (E n ) is bounded and has for upper-limit K, the inclusion E n C 
holds for n large enough. Hence cap(-E' n ) > cap(17) for every n. Therefore 

liminf cap(£' ra ) > cap(Q) . 
The open set VL being arbitrary, the desired conclusion holds. 
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QED 



Let fi be an open bounded subset of R , with S CC f2. We denote 
by .ffg(fi) the intersection sequence of the spaces i/g(fi n ) where (fi n ) is a 
decreasing sequence of open bounded sets, such that CC H n and £1 = 
n n Q n . One easily checks that Hq(£1) does not depend on the sequence (Cl n ). 

Lemma 3.3 Assume that \dd\ = 0. Then the following equality holds: 

cap(H) = inf J f \Vv\ 2 : v £ H^(U), v = 1 on S 
[Jn\s 

and there is a unique u £ Hq(Q) such that 

u = 1 on S and I |Vu| 2 = / |Vu| 2 = cap(fi) . 
Jr n \s Jn\s 

Moreover u is harmonic in tt\S and \ {u > 0}\O| = 0. 

Definition 3.4 Such a function u is called the capacity potential of with 
respect to S. 

Remark 3.2 

1. If dVt, is C ' , then the capacity potential u of Q with respect to S is the 
(classical) solution of (j5J) and is equal to the (classical) capacity potential of 
(see Remark 

2. In what follows, we study the energy of subsets ft DD S which is defined 
as the sum of the capacity and the volume of with respect to S (see ©). 
This energy is well-defined for bounded sets f2 DD S. It is why we assumed 
all the sets to be bounded. But let us mention that all classical results of 
this section hold replacing S bounded by Q\S bounded. We need this 
generalization in the proof of Lemma 14.51 

Proof of Lemma 13.31 The proof is easily obtained by approximation. By 
construction of cap(fi), we can find a decreasing sequence of open bounded 
sets Q n such that 

£1 CC Cl n , f] Q n = Q and cap(J7) = limcap(S7 n ) . 

n 

Let u n be the (classical) capacity potential of Q n . From the maximum 
principle, the sequence (u n ) is decreasing, and converges to some u which 
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is nonnegative with a support in 17 and equals 1 on S. In particular, {u > 
0} C 17 a.e. since \d£l\ = 0. Furthermore, by classical stability result, u is 
harmonic in 17 because so are the u n . Since we can find a smooth function 
4> with compact support in 17 such that = 1 on S, we have 

|Vu n | 2 < / |Vc/>| 2 = / |V</>| 2 , 
n n \s Jn n \s Jn\s 

which proves that (u n ) is bounded in H 1 (W N ). Thus the limit u belongs 
to H l (R N ). Since u n E H^(Q n ) with H^(Q n+1 ) C F^(n n ), « belongs to 
Hq (17 n ) for any n. Therefore u £ iTg(TT). In particular, the support of u lies 
in 17 = a.e.. So we have, 



cap® = limcaptO.Hlim/ \Vu n f 

[ |Vn„,| 2 > / \Vu\ 2 = [ \Vu\ 2 . (8) 
Jr n \s Jr n \s Jn\s 



= lim inf 

n 

For every n, 

cap(17 n ) = / |Vu n | 2 = inf \ / |Vw| 2 : u G (^n), « = 1 on S \ 
Jn n \s [Jn n \s J 

< inf \ f \Vv\ 2 : v G H^(T}), v = lonS\ , 
[Jn\s J 

since Hq(Q,) C ilgC^n)- Letting n go to infinity, we obtain 

cap(17) < inf J f \Vv\ 2 : v £ fl£(fi), u = 1 on 5 I . 

[Jn\s J 

From 0, we get the equality in the above inequality and the fact that u is 
optimal. Uniqueness of u comes from the strict convexity of the criterium. 

QED 

4 The discrete motions 

Let us fix h > which has to be understood as a time step. Let us recall 
that S is the closure of an open bounded subset of M. N with C 2 boundary. 
We introduce the functional space 

E(S) := {u e H 1 (R N )nL°°(R N ) : u = 1 on S and u has a compact support} 
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If S and S' are two compact subsets of M> N with C 2 boundary such that 
S C S', then we note that E(S') C E{S). 

For any bounded open subset O of with S CC O we define the 
functional : E(S) — > K by setting 

4(0, «) = / |Vu| 2 + 1 {U>0} f 1 + 

where is the signed distance to Q defined by ©, 1a denotes the indicator 
function of any set A C M N and r+ = r V for any r£R. We write J/i(f2, u) 
if there is no ambiguity on S. 

Let us recall some existence and regularity results given in [3]: 

Proposition 4.1 (Alt and Caffarelli |3j) Let Q, be an open subset ofR N 
such that Q\S is bounded and with Sccfl. Then there is at least a min- 
imizer u E E(S) to Jft,(0, •). Moreover u is Lipschitz continuous and is 
harmonic in {u > 0}\S. Finally 'H N ~ 1 {d\u > 0}) < +oo. 

Remark 4.1 We note that S CC {u > 0} because u is Lipschitz continuous 
with u = 1 in S. 

The existence of u and its Lipschitz continuity come from Theorem 1.3 
and Corollary 3.3 of 3 . The fact that u has a compact support is estab- 
lished in Lemma 2.8, and its harmonicity in Lemma 2.4. The finiteness of 
H N ~ 1 (d{u > 0}) is given in Theorem 4.5. 

We are now ready to define the discrete motions. 

Let f^o DD S be a fixed initial condition. We define by induction the 
sequence (f2^) of open bounded subsets of W N with J7 n DD S by setting 

Og := n and Q!^ +1 := {u n > 0} U {x € ft£ : danft^) > Mi 
where 

u n £ argmin jf v). 

vGE(S) 

We call discrete motion such a family of open sets. Of course, the dis- 
crete motion is defined in order that it converges to a solution of the front 
propagation problem Q (see Theorem 15.21 and Remark 14.2)1 . 

In order to investigate the behavior of discrete motions, we need some 
properties on the minimizers of Jh- 
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Lemma 4.2 Let f2 and u be as in Proposition \4- 1\ Let ft' = {u > 0} U 

where 



Q h :={y£n : d dn (y) >h}={yeR N : d s n (y) < -h} . (9) 
Then \dQ'\ = and u is the capacity potential to tt' 

Remark 4.2 We do not claim that u is positive in f2'. For instance, consider 
a set £1 with two connected components ill and such that S CC In 
this case, u = in £^2- Notice that it explains why we define := {ii n > 
0} U {x G 9^ : d dU h(x) > h}. Adding the set {x G Ojj : d dn h(x) > h} 
prevents the discrete motion from the sudden disappearance of a connected 
component. Indeed, the discrete motion is built in order to approach a 
solution of the front propagation problem (JTJ and a connected component 
which does not contain any part of the source is expected to move with a 
constant normal velocity —1. 

Proof of Lemma 14.21 Let us first notice that \dfl r \ = 0. Indeed we 
already know that \d{u > 0}| =0 (because its W^ -1 — measure is finite from 
Proposition 14. 1 I) . On the other hand dflh C {y 6 H : dont(y) = h} has also 
a finite Ti 1 ^^ 1 — measure thanks to 4, Lemma 2.4]. 

Let now e > be fixed and set, for any a > 0, Q a = {y G R N : do'(y) < 
a}. The set O a is open, bounded and satisfies Q' CC £l a . Moreover, since 
lf2 a - > 1q7 and O' is bounded with \d£l'\ = 0, for a > enough small, we 
have 

/n»,( 1 + K) + £e ' <10) 
Let v be the capacity potential of f2 Q and set 

v k (x) = v(x) + -d R N\p a (x) Vx G M. N . 

Then (v^) converges to v in H 1 (Wl N ) and \Q a \{vk > 0}| = 0. Therefore 



Jh{^,v k ) 



cap(O a ) + / ln Q (l + -c^) 



Since Jh(fl,Vk) > (fi, it), we get from (fTU|) 
cap(J7') > cap(f2 a ) 
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ln ( 1 + Uh 



> ]imJ h (n,v k ) - \ 

> J h (n, u )- ! i n „ (i + Uh) 

Jr n \s V h / + 

- L xs |v " |2 - in ^>»»( i+ H. 



> I \Vu\ 2 -e 



Thus J^n\s |Vn| 2 < cap(O'), which proves from Lemma 13.31 that u is the 
capacity potential of $7'. 

QED 

Next we need to compare solutions to Jh(Q, •) for different S and 0. 

Proposition 4.3 Lei S\ and S2 be the closure of two open bounded subsets 
of M> N with C 2 boundary, Q\ and Q2 be open bounded subsets of M. N such 
that S\ CC tt± and S2 CC Q2 - Let U\ and U2 be, respectively, minimizers of 
J, •) and Jfr 2 (Q,2, •)• ^/"^l C <S2 and f?i C O2, £/ten ui Ati2 and ui Vti2 
are, respectively, minimizers of J'^ 1 (Oi, •) and J^ 2 (^l2, •)• 

Remark 4.3 

1. In particular, if J^ 2 {Vt2,-) has a unique minimizer U2, then {ui > 0} C 
{u 2 > 0}. 

2. This Proposition still holds true if we replace, for £ = 1, 2, fij, bounded 
by Q.i\Si bounded; see Remark 13.21 and Lemma l4~5l 

Proof of Proposition 14.31 We have 

jf 1 (fix, tti Ana) + Jf 2 (^ 2 ,ni Vu 2 ) 
= J h 5l (^i^i) + 4 52 (^2,n 2 ) 

+ / (|V(«! A n 2 )| 2 - |Vm| 2 ) + (l {MlAU2> o } - 1{ UX >0}) (l + 
Jr n \Si ' V n J 

+ L» s (|v(ni v u2)|2 " |Vu2|2) + (i{uivu2>o} _ i ^>°} ) ( x + 
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Since Q\ C fi 2 we have < in R . Hence, a straightforward compu- 
tation leads to 



1 „ \ (. 1 



1{«iAm 2 >o} I 1 + T d h L ] + l{«iv« 2 >0} I 1 + v d 



< 1 



{«i>0} 



1 + + Wo} (l + ~^n 2 J 



Moreover, by classical results, 

|V(«i A u 2 )| 2 + |V(iti V u 2 )| 2 = |Vui| 2 + |Vu 2 | 2 a.e. in R^. 
It follows 

J^fn^U! Au 2 ) + jf 2 (0 2 ,"i Vn 2 ) 

< (S2 1 ,«i) + J? (n 2 ,«2) 

+ / (|V(«i A u 2 )| 2 - |V Ul | 2 ) + (l {uiAU2>0} - l{n 1>0 }) ( 1 + ^dfh 

But u\ A u 2 = ui on Si which gives 

J^(Q 1 ,u 1 Au 2 ) + J^(Q 2 ,u 1 Vu 2 ) < J^(n 1 ,u 1 ) + jf (n 2 ,u 2 ). (11) 

Since u\ and it 2 are minimizers we have 

Jfc x ( n i»«l) ^ jf^l.WlA^) and Jf 2 (0 2 ,u 2 ) < jf 2 (0 2 ,mVn 2 ). (12) 

The inequalities in (|11|) and l|12j) are therefore equalities. Hence ui Au 2 and 
u\ V u 2 are respectively minimizers of J h x (VL\, •) and J^ 2 (f2 2 , •). 

QED 

We define the energy £(Q) by 

£(H) = \U\ + cap(H). 

(compare with (jH). 

Lemma 4.4 Let (f^) fre a discrete motion with |<9f2g| = 0. Then the energy 
£(^n) is non increasing with respect to n. More precisely, 

£(n h n+1 )-sm < (i^\ { ^<_ M - i w >o } \ { ^<-m) l d k ^ °> 

where u n is a minimizer for ,7^(0^, •). 
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Proof of Lemma 14.41 Let us fix n. In order to simplify the notations, let 
us set 

Q := £1% , Cl h := {x £ 9 : d s Q (x) < -h} = {x £ Q : d m (x) > h} . 

Let uq be the capacity potential of 0, and u be a minimizer to J^(0,-). 
We finally set 9! := = {u > 0} U £l h . Recall that Q' G V and that 

1 90' | = 0: indeed this is true for n = from the assumption and by Lemma 
14.21 for n > 1. With these notations we have to prove that 



Then (u^) converges to uq in H (Mr) and {uj, > 0} = a.e. because 
{n > 0} C H and |<9ft| = 0. Hence 



£(«') < f (n) . 



For this we introduce for any k > 1 the function it^ defined by 




lim Jh(9, Uk) = lim 



k k 



R N \S 



Vn fe | 2 + l {ufe>0 }(l + rdri 




On the other hand, since cap(fi') = J*j 
|f2'| = |fJ'|, we also have 



Vu| 2 from Lemma 14.21 and since 




Writing that Jh(Q,,u) < Jh(9,Uk), we get the desired claim. 
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QED 

Next we show that the solution does not blow up when h becomes small. 

Lemma 4.5 Let R > and r G (0, R/2 1 ^ N ~^) be fixed. Let us also fix M 
such that + M > 4(iV — 2)/ro- Then there is some ho = ho(N,ro, R, M) 
such that, for any h G (0,/i ) and re (r , -R/2 1 /^- 2 )), /or any G open 
bounded, for any x G" wit/i r < dQ(x), i? < ds(x) and for any u minimizer 
to Jh(Q, ■), we have 

°Wo}un h ( x ) >r-Mh, 
where Qh is defined by (GJ). 

Proof of Lemma 14.51 The idea is to compare the solution with radial 
ones. For simplicity we assume that N > 3, the computation in the case 
N = 2 being similar. We also suppose without loss of generality that x = 0. 

Let us first investigate the problem of minimizing J h R (B£, •), where 
B r = B(0,r) and Br = B(0,R). Notice that neither the source B 1 ^, nor 
the subset B° is bounded but B^\B^ = Bji\B r is bounded so the previ- 
ous results on the minimization problem apply (see Remark IM. 2J) . Standard 




{u > 0} 



Figure 1: Illustration of the proof of Lemma 14.51 



symmetrization arguments show that a minimizer v to J h R (B^,-) must be 
radially symmetric. For p G (0, R), let us denote by v p the (radial) har- 
monic function which vanishes on dB p and is equal to 1 on 8Br. We also 



set Jh{p) '■= Jh R (Br,Vp)- Notice that a minimizer of J h R {B^,,-) has to be 
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either of the form v„ with p minimizer of Jh('), or constant equal to Vo := 1. 
Let us fix ho enough small in order that r + ft < R for ft G (0, fto). We have 

J h (0 )-J h (B r ,v )- hN{N + 1) , 

where oat_i is the volume of the unit sphere of WL N . For Jh{p) with p > 0, 
we distinguish two cases. If r + ft < p < R, then 

ajy_i(JV-2) 
~~ p 2 -^ - R 2 ~ N ' 

If < p < r + ft, then 

J h (p) 7V-2 l/(r + /») JV+1 p^ +1 (r + ft)/^ 



iV-2 , 1 /(r + Zi)^ 1 p w+1 (r + h)p N \ 
p 2-N _ R 2-N + I y AT(iV + i) + FT! iV J 



We show that vo cannot be a minimizer by comparing J/j(0 + ) with J/i(p) 
for < p < r + ft. Choosing p = with /? > 0, we have 

' -(J h (p)-J fc (0)) 



ajv-i 

N — 2 p 3 (r + ft)p 2 



JV-2 



+ 



l-(p/i?) iV - 2 ft(iV + l) ftiV 



jV-2 /3 3 ftV 2 r/? 2 \ 

" P \l_pN-2 h {N-2)!2/ R N-2 + # + 1 jV y 1 j 

Recalling that r G (0, i?/2 1 /( Ar " 2 )) is fixed, we choose 

> ggtg-q + i) 

and then ftp = ho(N, [3, ro, -R) > enough small such that 



1 - tt-^ > - and — < 1. 

i?^" 2 2 iV + 1 



(15) 



For all ft G (0, fto), we obtain that (|T3*|) is negative, which proves that vq is 
not a minimizer. 

Therefore minimizers have to be of the form v p for some p G (0, R). On 
(r + h, R), Jh(p) is increasing. For p G (0, r + h), we have 



4(p)_ (7V-2)V- W + p w (r + %/ 



a w _i (p 2 -^ - i? 2 "^) 2 ft ft 
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The stationary points of Jh on (0, r + h] satisfy 

(N-2) 2 1 
/(p) := \ * ±( r + h-p) = 0. (16) 

\p(l- (p/R) N ~ 2 )} 2 h y 

Notice that p \— > f(p) is convex on (0, r + h] and tends to +oo as p —* + 
and as p — > R~ . If we find some value p for which f{p) is negative, then 
there are exactly two solutions to (|16|) . 

For this, let us choose p = (3\/h with /3 > 0. Then 

Choosing /3 > satisfying ((11)) and 

< J 1/2 

and /iq satisfying (|15|) and 



,,,1/2 ^ r 

/?V < ( 17 ) 



we obtain that f(pVh) < for /i G (0, /io)- 

Let us fix /i 6 (0, /io) and let p\ and /32 be respectively the smallest 
and largest solutions to (|16[). With the arguments just developed above, 
we know that p\ < 0\fh < p2, where j3 is defined as above. Since J' h {p) = 
"Af-iP Ar_1 /(p) ) we have 

j'h(Pl) = UN-lPi' 1 f'{pi) 

-2(N-2) 2 (l-(N-l)( Pl /R) N - 2 ) 



N-l 



+ 



( Pl (l-( Pl /i?)^)) 3 
2(iV-2) 2 / . Ar . /pi\N-2\ 0- 



+ 



If we choose ho > satisfying Q15j). (|17j) and furthermore 

/ 1/2\ N ~ 2 

(»-l)f^-J <| and *j"<^. (18) 

we obtain that J' h \pi) < for h E (0, /io) and p\ is not a minimum to Jh- 
Therefore, Jh is increasing on (0, p\), decreasing on {p\, pz) and increasing 
on (p2, R). The minimum is achieved at p = p2- 
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Let us now estimate p2- We suppose that ho satisfies (|15|> . (fT7)l . (fT%)) and 

H) 16(iV-2) 2 
/io < 771T7 where 1 + M > v 7 . 

Then, for all /i G (0,/i ) andr G (r , i?/2 1/(7V ~ 2) ), we have r-Mh > r /2 > 
and we compute 



{N-2f 

(r - Mh) 2 {\ - {{r-Mh)/R) N - 2 f 



f(r~Mh) = - ' - (l + M) 



4(JV - 2) 2 . 

< o. 

Therefore p2 > r — Mh. 

To summerize, we know that, setting ho = ho(N, ro, R, M) small enough, 
for all h G (0, ho) and r G (ro, R/2 l / i - N ~ 2 ^), the problem consisting of min- 
lmizme Jh R {B c r ,-) has a unique solution w P2 , which is radially symmetric 
and such that p2 > r — M/i. 

Let now £ D, i ^ with i? < dg(x), r < du(x) and let u be a 
minimizer to Jft,(0, •). Since 5 C B c R (x) and $7 C Bf.(x), Proposition 14.31 
states that {n > 0} C {v P2 > 0} C B^_ Mh (x) (see Figure ^ for a picture). 
Since Qh C fl C B£, finally, we have di u> o}un h ( x ) — r ~~ 

QED 

Finally we explain that the set {n > 0} satisfies some inequalities in a 
viscosity sense. Here again the regularity results of Alt and Caffarelli [3] 
play a crucial role. Let £ be an open set with C 1 ' boundary such that 
S CC £ and X^S 1 is bounded. We denote by n^ the (classical) solution to 
(J3J (replacing ft by £), i.e., the capacity potential of £ with respect to 5. 

Lemma 4.6 Let Q, be a bounded open subset ofM. N with S CC SI and n oe 
a minimizer to J^(fi, ■). VFe se£ 

= {x G ft | dgn(x) > h} and ft' = {n > 0} U . 

Ze£ £ is an open bounded subset of~R N with C 1,1 boundary. 

1. [Outward estimate] Suppose that £ zs swe/t i/iai 

{n > 0} C £ and 3x G <9£ D <9{n > 0} . 
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Then 

1/2 



2. [Inward estimate] Let us now assume that S is sitc/i i/iai 
5CCS, Ecfl' and 3x £ <9£ n 90' . 

T/ten 



1/2 



Proof of Lemma 14.61 Let us set gn(x) = (1 + d s n (x)/h) + . We first prove 
the outward estimate. From |3J Lemma 4.10] we have 

u(x') / — -—- 
hmsup > y/gn{x) 

x' 6 {« > 0} 

for any ball B contained in {u = 0} and tangent to {u > 0} at x. Let 
v be the outward unit normal to S at x and r > be such that the ball 
-B := £>(x + ru, r) is tangent to S at x. Then 5 is also tangent to {u > 0} 
at x. Since by the maximum principle, u < Ug, we have 

it-, £ / \ i ,. u£(x') 

\vug{x)\ = hmsup 



x'—>x, x' '£ 



{«>o} rfs(x') 



u(x') 



> limsup 

as'-nc, x'e{«>0} a B\X ) 



We now turn to the proof of the inward estimate. We first prove that 
Ug < u in {ug > 0}. Indeed from Lemma 14.21 u is the capacity potential 
of Q'. In particular u is harmonic in Q'\S D u = Ug on OS 1 and 

= tig < u on <9{u|: > 0}. Hence < u in {it| > 0}. Let us note that 
u = on 90'. Therefore u{x) = ti^(x) = 0. 

We now consider two cases. If x ^ 9{u > 0}, then x G 90^; thus 
(1q(x) = — h and gn(x) = 0. But 0<M|<« = 0ina neighborhood of x so 
that Vu|;(x) = 0. Therefore 

\Vu§(x)\=0 = gn(x). 
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Let us now consider the case x 6 d{u > 0}. Then [3J Theorem 6.3] states 
that 

sup |Vu| < y/ gn{x) + m(r), 

B(x,r) 

where m(r) — > as r —* + . Since we want to prove that |Vug(x)| < 
y/gn(x), we can assume without loss of generality that Vn^x) / 0. Let 
v be the outward unit normal to S at x. Since v = — Vng(x)/|Vu^(x)|, 
for r > sufficiently small, the segment ]x,x — rv\ is contained in £ and 
in {ug > 0}, and thus in {u > 0}. So u is smooth at each point of this 
segment. Since moreover u > u$, we have, for some £ G (x,x — rz^), 

' u i( x — rzy ) — ii(x — rz;) = u(x) + (Vii(£), — ri/) 
< r (Vfi^Oc) + m(r)) . 



Therefore 



|V«i(x)| = r hm UsiX r TV) < y^(x) . 



QED 



5 Discrete motions and viscosity solutions 

Let us fix rio open and bounded such that S CC Slo- Let be a discrete 

motion with Qq = £Iq. 

Let us now introduce a lower and upper envelope for the sequences (f^) n 
as the time-step h tends to + : the upper envelope fC* is 

K*(t) := Let* : ^ ° + ' ^ 1 Xk £ ^' ) , (19) 
[ with Xfc — ► x and /ifcn/c — > i J 

while the lower envelope /C* is defined by its complementary: 

R N \K*(t) = LeR N : 3hk C ° + ' Hk Z t°°' ^ * ) • (20) 
[ with Xfc — > x and n^n/j — ► i J 

Lemma 5.1 The set /C* is closed while fC* is open. Moreover the maps 
t — > fC*(t) and t — > fC*(t) are left lower- semicontinuous on (0,+oo). 

Proof of Lemma 15.11 The fact that set /C* is closed comes from its 
construction since the upper limit of sets is always closed. The argument 
works in a symmetric way for /C*. 
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We now prove that t — > IC*(t) is left lower-semicontinuous on (0, +00) 
(see Section |21 for a definition). We proceed by contradiction assuming there 
exist t > 0, x G K.*(t), p > and a sequence t p — > t~ such that B(x,p) n 
/C*(t p ) = 0. Therefore (i c , (tp) (x) > p > for all p. 

Set i? = d s (x), r < min{p, i?/2 1 /( Ar ~ 2 )} and M > such that \/T+M > 
8(iV — 2)/vq. Then Lemma f4.5l states that there is some /io = ho(N, ro, R, M) 
with the following property: for any r G (ro/2, R/2 l l N ~ 2 ) and /i G (0, /io), 
for any SI with r < (1q(x) and for any u minimizer to J/j(0, •), we have 

d {«>o}uQ h ( x ) ^ r ~ M/l ' where ^ = { d f2 < ~ h }- 

For /i G (0, ho), let = [t p //i] be the integer part of t p /h. From the 
definition of K*(t p ) and ro, we can find some h\ G (0, ho) such that d^h (x) > 
ro for any h G (0, /ii). We are going to prove by induction that 

d Q h (x)>r - Mkh for all k G {0, . . . , A#}, (21) 

where /cq = [ro/(2M/t)j. Indeed inequality (|21(1 holds for = 0. Assume 
that it holds for some k < kft. Let u be a minimizer for Jh(^n h +kh^ ') anc ^ 
define 

n£ fc+ (*+l)ft = {n > 0} U {y G fi* h+fcA : ^ + Jy) > /i} . 

Then since ro — Mkh > ro/2 and ro — M/c/i < ro < R/2 1 ^ N ~ 2 \ we have 
from Lemma 14.51 recalled above that 

d Qh (x) > r - Mkh - Mh . 

n h +(fc+l)fc 

So (|2*Tj) is proved. 

Let us set r = ro/(4M) and fix s G (0, r). Let (fc/j) be such that k^h — > s 
as /i — > + . We notice that kh G {0, . . . , k^} for fa sufficiently small. Letting 
h — > + in inequality (|2*Tj) for any such (A; A ) implies that 

> r - Ms > r /2 > 0. (22) 

Since r does not depend on x and t p and since i p — ► i~, for p large 
enough, we have s = t — t p < r. Therefore, from (|22|).we obtain d^*^(x) = 
dfc*u +s)(x) > ro/2 > which is a contradiction with the assumption x G 
K*(t). 

The proof of the left lower semicontinuity of /C* is simpler. As above, 
we proceed by contradiction assuming that there exists x G JC*(t) for t > 
and a sequence i p — > t~ such that djp, Jx) > p > for all p. From the 
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definition of +1 , for (n/j) such that n^h — > i p and h sufficiently small, we 
have B p / 2 {x) C &n h - F rom the definition of £7^ +1 , we have therefore 

B p/2 „ h {x) C {y G Sl* h : d^^x) > h] C n£ h+1 . 
By induction we prove in a similar way that, for any k < p/(4/i), 

B p/2 - k h(x) c {y e fi£ h+(fc _i )h : d an^ +(Jk _ 1)k ( JC ) > /l l c • 
Letting now ft, — > + we get at the limit: 

Bp/4 0) n £*(tp + s) = for all s G [0, p/4] . 

Since p is independent of p, we get a contradiction by taking p big enough 
such that t — t p = s < p/4. 

QED 

Theorem 5.2 T/ie tu&e /C* (respectively /C*J is a viscosity subsolution (re- 
spectively supers olution) to the front propagation problem V = h(x, 0), where 

h(x,n) = -l + h(x,n) 

and h is defined by 

Proof of Theorem Let us set tl h := \J n {nh} x n n- Let (*o,a?o) G JC* 
with to > 0, be such that there is a smooth regular tube K r with K* C /C r 
and xq G dfC r (to). Without loss of generality we can assume that K,*C\dlC r = 
{(to,xo)}. Then by standard stability arguments (see [J]), one can find a 
sequence of smooth regular tubes K^. converging to JC r in the C 1,h sense (see 
Section |21 for a definition), and sequences hk — ► and — ► +oo such that 
S7 fefc C /C^, (nkhk,Xk) — > (io 5 £o)> £fc G <9fi^ and such that G dK.^(nkhk)- 
Let u be a minimizer to J^ J ,(f2^ fc _ 1 , •). By definition of the discrete mo- 
tion, we have 

n£j = {u > 0} U {y G : d« hfc (y) < -/»*} . (23) 

"fe- 1 

Let Ufc := u ^ r ^ nkhk ^ De the capacity potential of /C^n^/ifc). 

Let us first assume that Xk G > 0} for some subsequence of (xk) 
(still denoted by (xk))- The case Xk G int{u = 0} for any k is treated 
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later. From the discrete viscosity condition in Lemma 14, 61 and the inclusion 
^nt-i c fcr(( n k ~ l)^fc)j we know that 

|v« fc (a*)| > ^i + _^^( Xfe ) > (^i + _4, ((nfc _ l)/ifc) (x A; ) . 

Hence 

^4^((n fc -i)^)(^) < - 1 + \Vvk(x k )\ 2 . (24) 

Let us now recall that the normal velocity of /Cj? at a point (t, a;) 6 <9/Cjf 
is given by -J^ fc(t) (x). Since x fc G dK*(n k h k ), (n k h k ,x k ) -> (t ,so) and 

since /C£ converges to /C r , we have therefore that 

a 

4:*:(K-i)/i fe )( x fc) = d K*(n k h k )( x k) ~ h k^dl-k inkhk) (x k ) + h k e(k) 

where e(/c) — > as — > +oo and K^^) is the normal velocity of /C r at 
(to,xo). From (|24[1 we get for A; large enough, 

h{x k ,K k r {n k h k )) = -1 + |V^(x fc )| 2 > vg^j + e(fc) . 

Letting — ► +oo, we obtain 

h(x ,]C r (t )) = \\mh(x k ,K, k (n k h k )) > V^ XQ y 

The above equality is a straightforward application of |13| Theorem 8.33] 
since K k converges to K r in the C l,b sense (see Section |21 for a definition). 

We now assume that x k £ int{u = 0} for any A;. Then we have from (|23|) 
that 

Arguing as above we get 

-h k = d s hk (x k ) > d s m(nk _ 1)hk) (x k ) = h k V/£ } + h k e(k) , 
"ft- 1 

where e(k) — * 0. Dividing by h k and letting k — > +oo gives 

V£ X0) < -1 < -1 + |Vn^ (to) (x )| 2 = /i(xo,/C r (to)) • 
So we have finally proved that /C* is a subsolution. 
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We now show that /C* is a supersolution. The proof starts exactly as 
above: if there is a smooth regular tube K r with K r C /C* and some (to, Xq) £ 
<9/C* with to > and xo £ dJC r (to), then one can find a sequence of smooth 
regular tubes K% converging to K T in the C 1,h sense and sequences h k — > 
and — > +00 such that K^(nh k ) C for any n, (n k h k ,x k ) — > (to,xo), 
Xfe £ dVL hk (n k h k ) H dK^(n k h k ). Let u be a minimizer to ^(^t-D ')■ Then 
(ESI holds for 

Then using Lemma 14, 61 we get 



(25) 

where v k := u^ r ^ nkflk ' is the capacity potential of K^(n k h k ) with respect to 
S. Since G dVt hk (n k h k ), we have from (|23|) that d^ hk , h Jx k ) > — 

Therefore inequality (|25|1 can also be written as 

T-d K k {(nk _ 1)hk) (x k ) > -1 + \Vv k (x k )\ 2 . 
As before we have 

d k((n k -i)h k )( x k) = h kV^ xo) + h k e(k) . 

Hence 

h(x k ,K k r (n k h k )) = -1 + |V^(x fe )| 2 < vg^ + e(fc) - V*^ . 
Then we can complete the proof as above to get the required condition: 

h{xMt,))<v^ xo) . 

QED 

In particular we get immediately the following Theorem: 

Theorem 5.3 Let Oo be an open bounded subset ofM. N such that S CC Qq. 
Let K + and Kr be, respectively, the largest and smallest viscosity solutions 
to the front propagation problem with initial position Qq. Then 

Kr c K* c K* c K + . 

Ln particular, if the problem has a unique solution, i.e., K~ = K + , then 
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Proof of Theorem 15.31 Since JC + contains any subsolution and Kr is 
contained in any supersolution (see |S]), we have K,* C K. + and K.~ C /C*. 
Inclusion /C* C JC* holds by construction. Whence the result. 

QED 

6 The energy is decreasing along the flow 

Let fio be a bounded open subset of R N . Let we assume that the front 
propagation problem (""J with initial position Qq has a unique solution and, 
furthermore that 

\dn \ = and |/C + \/C~|=0, (26) 

where KL + and K~ denote the maximal and minimal solutions respectively. 

Theorem 6.1 Under assumption \26]) . there is a set T C [0, +oo) of full 
measure such that 

£ (/C+(t)) < £ (/C+(s)) for all s,t€T, s<t . 

Remark 6.1 Assumption (|26|) is not too restrictive. Indeed, it is generic in 
the following sense: let (J7q)^ > q be a strictly increasing family of bounded 
open initial positions containing the source, i.e., 

for all < A < A', S CC fi£ CC . 

If JCt (respectively JCT) is the maximal (respectively minimal) viscosity 
solution to (P) with initial position Qq, then ((2T|) holds for all A > except 
for a countable subset. See |H] for details. For simplicity of notations, we 
have chosen to consider the case A = 1 and to assume that ()26|) holds for 
the initial position 0,q. 

Proof of Theorem 16.11 Let (^) be a discrete motion starting from Qq. 
Recall for later use that, from Lemma 14.41 

S(pk) < £(Th) Vn > 0, \fh > , (27) 

because we have assumed that \dQo\ = 0. Let K,* and /C* be the associated 
generalized evolutions defined by (|T^|) and 1(2*0)) , We have 

KT c /C* c K? c JC + . 
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Let 

T := {t G [0, +00) : |/C+(t)\/C-(t)| =0} . 

From assumption (|26|) and Fubini Theorem, the set T is of full measure in 
[0,+oo). 

We first prove that 

S(K+(t))<S(n^ VtGT. (28) 

For this, let t G T, /i^ — > + and — > +00 such that h^n\~ — ► 
For simplicity we set := Since the Kuratowski upper limit of the 

(f2fc) is contained in /C + (t), which is a compact subset, the sequence (O^) 
is bounded. Since moreover the upper limit of (R^\fifc) is contained in 
M N \JC~ (t), the latter with a boundary at a positive distance from S, there 
is some r > such that 5V C for any k sufficiently large (see (J7J) for a 
definition of S r ). Since finally the capacity is non increasing with respect to 
the inclusion, we get from Lemma 13.21 

liminf cap (T^T) > cap (£+(£)) . (29) 

k 

The next step towards (|28j) amounts to show that 

\IC + (t)\ < liminf \n k \ . (30) 

Let R > be sufficiently large so that fC + (t) CC B R , where B R = B(0, R). 
By definition of the Kuratowski upper limit and the construction of /C* , we 
have 

lB R \K,(t) > lim supl Bfl \n fc . 

k 

Fatou Lemma then states that 

\B R \lC{t)\ > Iimsup|Sfl\n fc | , 

whence (|3H)l since K.~(t) C /C*(i) and |/C + (i)| = |/C~(i)| because t G T. 
Combining ((221), © and (J23) finally gives 

< liminf £(fi^) < £(Th) Vi G T . 

A; 

This proves (|3H|). 

Let now < s < t with s,i G T. From the uniqueness of the solution 
starting from Kq, the maximal solution to the front propagation problem 
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starting at time s from IC + (s) is equal at time t to fC + (t). Since \dlC + (s)\ = 0, 
because s £ T, inequality ()28j) states that 

£ (JC+(t)) < £ (£+(*)) , 

which is the desired result. 

QED 
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